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We study states in the superconducting gap induced by magnetic impurities using self-consistent 
quantum Monte Carlo with maximum entropy and formally exact analytic continuation methods. 
The magnetic impurity susceptibility has different characteristics for To ^ Tco and To > Tco (Tq: 
Kondo temperature, Tco: superconducting transition temperature) due to the crossover between a 
doublet and a singlet ground state. We systematically study the location and the weight of the gap 
states and the gap parameter as a function of Tq/Tco and the concentration of the impurities. 
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Introduction. The problem of magnetic impurities in 
superconductors has been of great interest for a long 
time. It involves the competition of two distinct physical 
phenomena. The superconducting state, characterized 
by a transition temperature Tco, is formed by the coher- 
ent pairing of electrons with time-reversal symmetry. On 
the other hand, the magnetic moment on the impurity, 
characterized by a Kondo temperature Tq, couples anti- 
ferromagnetically to the conduction electrons and breaks 
the time-reversal symmetry of the conduction electrons. 
Thus even a small amount of magnetic impurities can 
break the pairs and form states within the superconduct- 
ing gap. 

Although many theoretical attempts have been 

made to describe these gap states induced by mag- 
netic impurities in superconductors, most of them are 
based upon perturbative approaches specialized for ei- 
ther To <^C Tco or To 3> Tcq. For example, Zittartz, 
Bringer, and Miiller-Hartmann (ZBM) [§ studied the 
problem using the Nagaoka-Suhl approximation which 
is a self-consistent high-temperature perturbation theory 
in the exchange interaction J which becomes exact when 
To <IC Tco- Matsuura |^ used an approximate interpo- 
lation method between two regimes, and found that for 
a single spin-i impurity system the ground state (gap 
state) is a singlet (doublet) for Tq ^ Tco while the ground 
state (gap state) is a doublet (singlet) for To -C Tco. It 
is easy to understand why the impurity has such states 
in these limits. When To is vanishingly small (J = 0) 
the impurity decouples from the superconducting host 
and has a doublet ground state. Even when Tq is finite, 
but small (To <C Tco), the superconductivity forms before 
the impurity can be screened by the Kondo effect, so it 
maintains its doublet ground state. However, in the limit 
where To ^ Tco, the Kondo effect completely screens the 
impurity before the superconducting transition temper- 
ature is reached. Thus, the impurity will have a singlet 
ground state in this limit. 

Nonperturbative methods have also been employed. 
Satori et al. applied the Numerical Renormalization 
Group (NRG) method to the problem of a single spin-i 
impurity in a BCS superconductor at zero temperature. 
This, however, does not account for the self-consistent 



effect of a finite concentration of impurities or address 
thermodynamic properties. Jarrell et al. introduced 
a self-consistent Quantum Monte Carlo (QMC) method 
to the dilute magnetic alloy superconductors and used 
the Maximum Entropy Method (MEM) ||| to obtain the 
superconducting density of states. Despite the fact that 
this QMC method treats the impurity nonperturbatively, 
the MEM alone does not produce enough resolution for 
a systematic study due to the complicated structure of 
the density of states. 

In this Letter we discuss a novel method to over- 
come these limitations. We use a combination of a 
self-consistent QMC method for the magnetic impu- 
rities, Eliashberg equations for the host superconduc- 
tors, MEM for the analytic continuation of the impurity 
Green's function, and a formally exact analytic continu- 
ation method (?) for the host Green's function. 

Formalism. We model the magnetic impurity with a 
symmetric Anderson model in the limit of infinite metal- 
lic bandwidth. This model is characterized by an on- 
site Coulomb repulsion U and a hybridization width 
r = uNqV'^, where V is the hybridization matrix ele- 
ment and A'o is the density of states at the Fermi energy. 
In the limit U ^ T , a, spin-^ magnetic moment on the 
impurity couples antiferromagnetically to the conduction 
electrons with an exchange J = — 8r/7riVoJ7. Then the 
Kondo temperature Tq is defined by jsj 
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and it is related to Tk defined in high-temperature per- 
turbation theory in J by Tk/Tq ~ {2/TTf 

We describe the host superconductor with the self- 
consistent Eliashberg equations. We assume that the 
conduction electrons interact with phonons with a cou- 
pling strength A(0) and the phonon density of states T(ijj) 
consists of a Lorentzian-like peak centered at wo- With 
the normalization Jg°° F{u))du! = 1, the electron-phonon 
spectrum function, a'^F{uj), is given by 
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When the width lol = 0, it represents Einstein phonon 
with frequency ojq. In order to have numerical stability 
and a continuous density of states, we choose a finite 
value oi ijJl- 

We assume that the concentration c of impurities em- 
bedded in the superconductors is small, and there are no 
correlations between them. In this dilute limit, we can 
describe the effect of the impurities on the superconduct- 
ing state using a standard perturbation theory with an 
average t-matrix approximation. 

We use 2x2 matrix Nambu formalism to represent 
both the host and impurity Green's functions. The QMC 
simulation determines the fully-dressed impurity Green's 
function starting from the initial, [/ = 0, impu- 
rity Green's function G°. Then we renormalize the host 
Green's function with Gd, and update G°. We iterate 
this self-consistent set of equations ||] until convergence 
is reached. 

This algorithm produces the fully-dressed impurity 
Green's function G(i(T) in imaginary time as well as the 
gap function A(ztt'„) and the renormalization function 
Z{iLo) in Matsubara frequencies u;„ = (2n -I- 1)'kT . Thus 
proper analytic continuations to real frequencies lj are 
required to obtain the density of states and the impu- 
rity spectral function. Since it has been shown that the 
MEM 1^ gives good results for various impurity prob- 
lems, we use the MEM for the impurity Green's func- 
tion to obtain the impurity spectral function A(cj) and 
Gd(uj). It was necessary,however, to develop a technique 
to analytically continue the non-positive definite off- 
diagonal term (A([j))i2, which we will discuss in a sub- 
sequent paper. Marsiglio et al. 0| presented a formally 
exact analytic continuation for the Eliashberg equations. 
We extend their formula to this problem by including 
the diagonal part (G(i(w))ii and the off-diagonal part 
(Gd(w))i2 of the impurity Green's function, and we have 
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Properly taking care of the singularities in Eqs. and 
® 101' solve these self-consistent equations with 
A(i[x)„) and Gd{uj). Then the superconducting density 
of states N{u!)/Nq is given by 



Re 
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Results. We have systematically studied the prob- 
lem for many different set of parameters. We choose 
Tco — 0.20 particularly for the results in Figs. such 
that LOa = 0.75, A(0) = 3.5 and ujl = 0.25. First of 
all, we study a single impurity in a superconductor. We 
may infer the form of the magnetic impurity suscepti- 
bility Ximp for small T based on Matsuura's observation 
regarding the crossover between a doublet and a singlet 
ground state, and we have 



T for a doublet ground state 
Te*^/-^ for a singlet ground state 
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where S is the energy difference between the ground state 
and the first excited state (gap state). Our QMC results 
for the inverse magnetic impurity susceptibility Xj~ip(T) 
are shown in Fig. ^ for various values of Tq/T^q- With 
extrapolation to T — > 0, X^p '^^^^ intercepts when 
Tq/Tcq ^ 1 (the open symbols) but not when Tq/Tco ^ 1 
(the filled symbols). With Eq. (|^), this implies that the 
doublet ground state is favorable for Tq/Tco ^ 1 and 
singlet for To/Tco ?J 1. It should be possible to see this 
behavior in NMR measurements of I/T2 oc Ximp- 

When the concentration of impurities is small but fi- 
nite, their effects can be seen in the dynamics of the 
problem. We use the rescaled concentration of impuri- 
ties c* = c/[(27r)2A'^gTco] and choose T = 0.05 where the 
superconducting gap of the pure host is 95% of its maxi- 
mum value. Fig. ^ shows the superconducting density of 
states for various values of Tq/Tcq when c* = 0.03. For 
energies above the gap the phonon induced structure in 
the density of states is the same as for the pure super- 
conductor (dotted line) but there is also structure in the 
gap. As To increases through Tco in Fig. ||, the gap states 
move toward the center of the gap and move back to the 
gap edge. The gap states form at the center of the gap 
at Tq/Tco — 1 where the crossover between doublet and 
singlet ground state occurs and the pair-breaking rate is 
largest due to the mixing of these states. Fig. ^ shows the 
position of the gap states (relative to the center of the 
gap) for various values of Tq/Tco and the corresponding 
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ZBM results. There is a much better fit if we use Tq in- 
stead of Tk in their formula. This is, in fact, the param- 
eterization chosen when their formula was fit to different 
experiments Even so, for Tq/Tcq < 1, the gap states 
move faster than ZBM's as a function of Tq/Tcq, which 
also agrees with experiments . 

We find that the weight of the gap states (not shown 
here) is linearly proportional to c*, which is consistent 
with the experiments and the ZBM result. However, 
while the location of the gap states is a universal function 
of Tq/Tco, the integrated weight also depends upon the 
couphng A(0) (increasing with decreasing A (0)). Further- 
more, for fixed A(0), the spectral weight of the gap states 
as a function of Tq/Tcq is not symmetric with respect 
to log'^jTo/Tco) which is in agreement with the NRG re- 
sults H). When Tq/Tco ^ 1 and the impurity ground 
state is a doublet, we find that the weight is relatively 
flat. On the other hand, when Tq/Tcq }^ 1 and the impu- 
rity ground state is a singlet, the weight increases with 
To/Tco- 

To understand the asymmetry note from Eq. (|^) that 
any contribution to ImA(ti;) at frequencies lower than 
the gap, will yield states in the gap, whether this con- 
tribution is from the impurity (i.e. —cV'^{Gd{(^))i2), or 
the remaining host portion of Eq. (^. At low frequen- 
cies when Tq/Tco < 1 both the real and imaginary parts 
of the off-diagonal impurity Green's function have oppo- 
site sign compared to the host result; i.e., the real and 
imaginary parts of the impurity gap function are out of 
phase with the host result. Hence the impurity and host 
contributions to the gap states tend to cancel. When 
Tq/Tco > 1, they have the same sign and tend to add. 
Thus, when the impurity ground state crosses over from 
a doublet to a singlet, it is accompanied by a phase shift 
of the impurity order parameter ||l2|| and an asymmetry 
in the weight of the gap states. 

As c* increases, the superconducting gap parameter Aq 
and the transition temperature Tc for the magnetic alloy 
superconductors are reduced from the values for the pure 
superconductor Aqq and Tco respectively. Fig. ^ shows 
that for small c* the initial reduction of Aq is largest 
when To/Tco ~ 1. In Fig. ||, we plot Aq/Aqo vs. c* for 
different values of Tq/Tcq. For larger c* the reduction 
of Aq is again largest for Tq/Tc ~ 1, where Tc is the 
transition temperature at that value of c*. The steepest 
slope of Aq/Aqo is obtained at larger concentration of 
impurities for Tq/Tcq ^ 1 and at smaller concentration 
of impurities for Tq/Tcq ^ 1. This may be understood 
as follows: When Tq/Tcq < 1, increasing c* decreases Tc 
and thus increases Tq/Tc toward 1, and hence increases 
the pair breaking rate. Therefore, when Tq/Tcq < 1, the 
curvature of Aq/Aqq is initially negative. In the other 
regime, Tq/Tcq > 1, the initial curvature of Aq/Aqq is 
positive. 

Note that the maximum initial slope |(9Aq/9c)^^q| 
occurs when Tq ~ Tcq, but the maximum initial slope 



|(9Tc/9c)^^q| occurs when Tr- ~ Tco These facts 

are independent of the coupling strength A(0). This im- 
plies that Aq and Tc have different functional forms even 
to the first order in c; i.e., Aq/Tc is not constant. 

Conclusions. Using a self-consistent QMC, MEM and 
a formally exact analytic continuation method 0, we 
have calculated the density of states for a dilute mag- 
netic alloy superconductor. We notice that the mag- 
netic impurity susceptibility has different characteristics 
for To/Tco ^ 1 and Tq/Tco ^ 1 due to the crossover 
between a doublet and a singlet impurity ground state. 
The gap states form at the center of the gap when 
Tq/Tco (not Tk/Tco) ~ 1 where the crossover between 
the doublet and singlet ground state occurs. Accom- 
panying this crossover there is a phase shift of the im- 
purity order parameter, and the spectral weight of the 
gap states as a function of Tq /Tco is not symmetric with 
respect to log^(ro/rcQ). As the concentration of impu- 
rities increases, the superconducting gap parameter Aq 
is reduced with different curvatures for Tq/Tco < 1 and 
Tq/Tcq > 1. We also find that Aq /Tc for fixed A(0) is 
not constant but a function of both the concentration of 
impurities and Tq/Tcq for a dilute magnetic alloy super- 
conductor. 
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FIG. 1. Inverse magnetic impurity susceptibility xTmp ''S- 
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the magnetic moment /i^ and To. 



T/To for various values of Tq/Tcq. Ximp is normalized with 
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FIG. 2. 
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Superconducting density of states N{uj)/No for 



various values of To/T^o when c* = c/[(27r)^iVoTco] = 0.03. 
Since N(uj) is symmetric, only positive frequencies arc plot- 
ted. The dotted line is for the pure superconductor with 
Tco = 0.20 and T = 0.05. Aqo is the superconducting gap 
parameter for the pure superconductor. 
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FIG. 3. Normalized location of gap states, y = <5/Ao, for 
various values of To/Tco when c* = c/[(27r)^AroTco] = 0.01. 
6 is the energy difference between the center of the gap and 
the peak of the gap states, and Ao is the superconducting 
gap parameter. The filled circles are our QMC results. The 
dashed line is the ZBM result with Tk and the solid line with 
To. Multiple points appear when y = since when To ~ Tco 
the distributions of gap states for positive and negative ui 
combine into a single peak centered at a; = 0. 
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FIG. 4. Ao/Aoo vs. c* for different values of Ta/Tco- Ao 
is the superconducting gap parameter for the magnetic alloy 
superconductors and Aoo is for the pure superconductor when 
T = 0.05. 
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